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Abstract 

Mh J We investigate the collective excitations of a single-species Bose gas at T = 

in a harmonic trap where the confinement undergoes some splitting along one 
spatial direction. We mostly consider onedimensional potentials consisting of 
two harmonic wells separated a distance 2zq, since they essentially contain all 

f- >, , the barrier effects that one may visualize in the 3D situation. We find, within 

a hydrodynamic approximation, that regardless the dimensionality of the sys- 
tem, pairs of levels in the excitation spectrum, corresponding to neighbouring 

(^ ■ even and odd excitations, merge together as one increases the barrier height 

ov 

On 



up to the current value of the chemical potential. The excitation spectra 
computed in the hydro dynamical or Thomas-Fermi limit are compared with 
the results of exactly solving the time-dependent Gross-Pitaevskii equation. 
We analyze as well the characteristics of the spatial pattern of excitations of 
threedimensional boson systems according to the amount of splitting of the 
p* condensate. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein condensation (BEC) |l|-§[ has inspired 
an impressive amount of further laboratory research together with important theoretical 
advances. The overwhelming material produced in the last three years has been recently 
reviewed by several authors @-§|. Among the various fascinating subjects that have at- 
tracted the attention of the community, the predictions of the spectrum of low-lying ex- 
citations and sound velocity |7||| in a hydrodynamic picture based on perturbations over 
the Thomas-Fermi regime, and their experimental verification PHT^fl play a special role, 
since the excellent agreement between theoretical and observed figures bring into evidence 
the validity of the mean field description of large condensates. Methods to obtain exact 
solutions of the time-dependent Gross-Pitaevskii equation (GPE) have been developed |13 



which permit further confirmation of the accuracy of the TF+hydrodynamic approach in the 
thermodynamic limit; moreover, the large amplitude motion of a condensate has been in- 



vestigated in a fashion that demonstrates the validity of a scaling hypothesis [14], according 



to which a condensate may evolve mantaining the appearance of a barometric distribution 
in a quadratic field. 

A particular subject of the BEC field is the possibility of forming and studying double 
condensates from several viewpoints. Starting from the original MIT trap [§] and its later 
improvements that allowed, for instance, the observation of interference fringes and the 
experimental realization of the coherence of matter waves |L5], the isomorfism between 



atomic condensates and optical devices has been intensively investigated |nj. Tunneling 
between the wells and the features of Josephson currents now constitute common subjects 
in the literature ]T7| , |T%| |. A special chapter includes the binary mixtures of condensates or 
double condensates of distinguishable particles, typically associated to hyperfine splitting 
of the magnetic atom levels in the trap jl9| . In particular, a mean field description of 



two species, spatially split double condensates has been advanced in Ref. J23[. Indeed, 
most theoretical approaches regarding the low energy dynamics of a double well pO| , ^T| , 
tunneling and/or Josephson effects consider essentially this type of distinguishable systems 
- or concentrate upon dynamical aspects that exclude the quantum and thermodynamic 
analysis of the double well as a subject in its own right. The investigation of the properties 
of single-species, double condensates with substantial spatial overlapping could, however, 
become relevant in the near future in view of the current possibility of adiabatically distorting 
traps |22|] creating long-lived systems with various minima and barriers [|32 |. 



Recently, we have advanced some results on the thermodynamics of BEC of indistiguish- 



able particles in harmonic double wells |2J], which indicate that with a minimum amount 
of complication with respect to the well-known harmonic traps, one can trace the evolu- 
tion of the atomic spectrum as a single potential becomes progressively split into two with 
the same frequency, at a distance 2z apart along the z-axis, and separated by a barrier 
Vo = Mu^Zq/2. The major features of BEC, under such conditions, can be accurately 
described and interpreted in terms of a smooth transition from a system with N particles 
to two identical ones with N/2 each. With this in mind, in the present work we intend 
to analyze the effects of the presence of such barriers upon the collective excitations of the 
condensate. For this sake, in Sec. 2 we concentrate on the onedimensional, harmonic double 
well problem; we extract the exact eigenspectrum of the linearized hydrodynamical equa- 
tion |7| , and confront the results thus obtained with the exact Fourier spectrum of the GPE 
solutions. In Sec. 3, we examine the threedimensional system and examine, in particular, 
the spatial patterns of the perturbation as the well splitting increases. Section 4 is devoted 
to the full discussion of the results here presented. The final summary and comments are 
the subject of Sec. 5. 

II. ONEDIMENSIONAL BOSE SYSTEM IN A HARMONIC DOUBLE WELL 

The general formalism leading from the time dependent Gross-Pitaevskii equation (GPE) 
to the description of the low energy excitations of a Bose-Einstein condensate has been 
derived by Stringari |F] and reviewed, for example, in Refs. [f|||. The procedure consists 
of linearizing the hydrodynamical version of the GPE with respects to small deviations Sp 
from the Thomas-Fermi equilibrium density inside the condensate 

PTF(r) = ^kM (2.D 
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where p is the chemical potential evaluated in the Thomas-Fermi regime and g = Airh 2 a/M 
is the hard-sphere interaction strength for an s-wave scattering length a. This linearized 
equation reads 



MQ 2 5p = -V.[(p~V ext )V5p] 
Its one dimensional version, namely 



Mn 2 5p = —^- 

oz 



d 
(fJ'-V ext (z)) — Sp 



(2.2) 



(2.3) 



with V ext a harmonic potential with frequency u, after an adequate scaling is precisely 
the differential equation satisfied by the Legendre polynomials P n (JMuj 2 /2poz), the cor- 
responding eigenvalues being of the form Q^ = u 2 n(n + l)/2 pB| . Furthermore, the TF 
chemical potential po reads, in this case, 



p = — - 6vr N — 
2 V a. 



(2.4) 



with a z = Jh/Mu the corresponding oscillator length. 

Our current purpose is to find the exact solutions of Eq. (|2.3| ) when V ext is the double 
well potential 



V ext = -Muj 2 {\z\-z ) 2 



(2.5) 



In this case, within the TF approximation, the ratio of the chemical potential p to that of 
a single well is found to satisfy the relationships 



/' 



Ho 2 2/3 



if V >p 



(2.6) 



3/2 + 3 (vX 12 fpT 
p J 4 V Mo/ VAk), 



2 \m 



3/2 



1 otherwise 



where the barrier effect is fully included in the ratio Vo/po- The decrease p/po of the chemical 



potential is plotted as a function of the dimensionless variable JV /p = z /Z , with Z the 
single-condensate size, in Fig. 1, while in Fig. 2 we display three TF densities, respectively 
corresponding to p = 5/2 Vq, Vq and 1/4 Vq. 



We now look for the solutions of Eq. ( |2.3| ) in a fashion that closely parallels the com- 



putation of the eigenvalues of Schrodinger equation for the double well hamiltonian |25|] . In 
other words, we look for eigenvalues Q 2 , and eigenvectors Sp u of the form 



u> 



ni = —v(v + \) 



(2.7) 



and 



5pp 



Pu 



yjMu*/2y.(\z\-zv) 



(2.8) 



with P v (x) a generalized Legendre function |27[. To obtain the complete solution of the 
hydrodynamic problem we observe that i) the operator on the right-hand side of ( |2.3|) has 
well defined parity and ii) the boundary condition that guarantees real values of Q is the 
nondivergence of the density fluctuation on the edges of the condensate [p8 |. Taking this 
into account, we construct the eigensolutions ( |2.8[ ) separately for positive and negative z 
and request the matching conditions at z = [p5[ . 



Pu 



-^Muu 2 /2pz 



for odd solutions, 



(2.9) 



Pi. 



-JMu 2 /2pz 



for even solutions. 



(2.10) 



When the barrier height reaches the current chemical potential value p, Eq. ( |2.3| ) is well 
defined in two separate domains (cf. Fig. 2) and we can construct the density fluctuations 
as even and odd superpositions of the exact single-well solutions. 



(fyeven = A/" {Sp (z - Z ) Q(z) + 5 p Q (- Z - Z ) Q(-z)} 

Spodd = A/" {Sp (z - z Q ) Q(z) - Sp (-z - z ) Q(-z)} 



(2.11) 



being 5po one eigensolution for the ID harmonic trap with a length scale corresponding to 
halving the particle number, and B(x) the Heavyside function. Thus, once the condensate 
is split into two pieces, under the TF approximation, the excitation spectrum for any finite 
barrier Vq higher than the ground state energy p is the single well one, however doubly 
degenerate. The appearance of the ID excitation spectrum as a function of the separation 
Vq/ p = zq/Z is depicted in Fig. 3; moreover, to inspect the departures in the behavior 



of the excitations from the hydrodynamic limit, we have integrated the time dependent GP 
equation and display with symbols some indicative results. This integration is achieved 
in two steps: we first build up the condensate wavefunction solving the time-independent 
problem by the steepest descent method |[30|j , and afterwards we add some noise to the 
condensate density po- We then follow the time evolution of the perturbed density and 
compute the fast Fourier transform (fit) at one minimum of the double well potential, which 
provides the spectrum of small amplitude oscillations of the double condensate H31f| . This is 
illustrated in Fig. 4, where we plot this spectrum for four choices of the ratio JV /p and 
an interaction strength Na/a z = 10. We realize that for very low barriers, the spectrum 
is the one corresponding to a single well, but as one increases the separation between the 
minima, two quasidegenerate branches can be distinguished. These branches completely 
merge at infinitely large barriers, however, they remain observable at finite values of Vq due 
to the limited resolution on the fit calculation. We notice that although the slopes of the 
odd eigenvalues of the TF hydrodynamics exhibit a sharp discontinuity at the barrier limit 
Vo/p = 1, the fft results behave smoothly across the transition. We shall analyze this issue 
later in Sec. 4. 



III. THE THREEDIMENSIONAL BOSE SYSTEM 



We now discuss the solution of the full equation (|2.2| ) for a cylindrical potential of the 
form 

Tr , . MoJi 9 Moot ,. n9 . 

V ext (r) = —^L r 2 + —^L {\z\ - z ) 2 (3.1) 

with a given aspect ratio A = uj z /uo r . In this case, the chemical potential \x in the Thomas- 
Fermi limit is related to the single-well one given by 

Ao = -(l5A^) (3.2) 

with Cj = (uj 2 uj z )3 and a the corresponding oscillator length, according to 
H 1 



if V > jj, 




The ratio /x//xo is displayed in Fig. 5 as a function of the dimensionless variable JV Q /n Q 
and in Fig. 6, we display three contour plots of a TF density - i.e., the condensate edge 
and two inner equiprobability contours - for a number of particles such that pi = 5/2 Vq. 
Alternatively, these curves represent three different condensate edges for /j, = 5/2 V , V and 
1/4 Vq, in correspondence with the ID density profiles of Fig. 2. 

Due to the different definition of the z— components of the restoring force for positive 
and negative values of the coordinate, standard methods p|j8| j28|j2m employed to derive 
the eigenfunctions and eigenvalues for anisotropic single-well potentials are not especially 
useful in this case. We have rather found convenient to expand the solutions 5p(r, z) in the 
spherical basis p( 2n >{r) r l Y[ m (9, 0)) 0, and diagonalize the matrix representing the operator 
V. [(V ext — /i) V] jM. The eigenspectra thus obtained are labelled by the magnetic quantum 
number m and in order to illustrate the appearance of the spectra Q m k, with k just numbering 
the level for the given m, in Fig. 7 we display these energies for (a) m = and (b) m = 1 
as functions of the well separation. We appreciate the general trend similar to that already 
analyzed in the ID case, namely, the gradual merge of pairs of levels corresponding to even 
and odd eigenfunctions and their collapse into the doubly degenerate spectrum of single- 
well excitations. The apparent discontinuity of the odd level energies when the chemical 
potential equals the barrier height is also visible; furthermore, the eigenfrequencies arising 
from the numerical diagonalization can be compared with those provided by the Fourier 
transformation of an initial perturbation that evolves according to the full GPE dynamics. 



Figure 8 displays a typical Fourier spectrum for a weak barrier with JVo/fX = 0.14; some 
eigenvalues of the dynamical matrix are indicated on the abscissa axis so as to remark 
the good agreement between both procedures, in spite of the inaccuracies intrinsic to the 
Thomas-Fermi approximation. However, the way in which the excitation was created in 
the fit method - i.e, adding some random fluctuation to the TF ground state - does not 



allow us to control exactly the particular selection of the fluctuating modes or their relative 
intensity. Thus, there exist eigenfrequencies of the problem that we cannot clearly see in 
the fit spectrum of an arbitrary initial condition like the typical example in Fig. 8. 

The shape of the perturbation 5p m k(r,z) as a function of the coordinates inside the 
condensate deserves some comments. We have chosen to analyze the transition density 
patterns in the following way: keeping Fig. 6 as a reference, we consider, for a given 
condensate characterized by a choice of the chemical potential p, the equiprobability contours 
Ptf{^, z) = Po with pq a constant that vanishes at the condensate edge and increases towards 
the center of the trap. In Fig. 9, we are drawing four portraits corresponding to the first odd 
parity excitation for m—0, for different well separations (or to different chemical potentials 
for fixed barrier height). Each of these sets contains curves Sp[r(z),z] as functions of z for 
a selected equiprobability contour p . Full lines indicate the condensate edge po= , while 
the different dashed lines correspond to inner contours of increasing equal density. Figures 
10 to 12 respectively display the same patterns for the first even parity excitation with m=0 
and for the first odd and first even parity eigenmodes with m=l. Very similar scenarios 
are obtained for other low energy levels, the present choice being convenient to observe 
the evolution of the perturbation and its nodal characteristics in terms of relatively simple, 
however rich enough, spatial structure. We examine these results further in Sec. 4. 

IV. DISCUSSION OF THE RESULTS 

From the previous chapters, we recognize that the most important differences between 
the present collective dynamics and the well-known excitation spectra of Bose condensates 
in anisotropic harmonic wells, can be fully understood in terms of the dynamical properties 
of the ID double well. It is interesting to remark that, in this situation, the hydrodynamical 
wave equation ( |2.3|) can be solved almost analytically in the same spirit, and with the same 
degree of complexity, as the quantum problem of the Schrodinger eigenspectrum of a har- 
monic double well. The behavior of the excitation energies under a progressive well splitting 
(cf. Figs. 3 and 6) is a monotonic decrease together with an evolution towards doubly 
degeneracy of the excitations of a single harmonic well. This can be understood in terms of 
a simple classical picture that illustrates the effect of the barrier; in fact, the hydrodynamic 
equation (|2.3|) may be interpreted as the classical wave equation of a perturbing density fluc- 



tuation 5p travelling inside the condensate with a local velocity v(z) = J(p — V ext (z)) /M. 
With this in mind, we can estimate the time r employed by the perturbation to go from one 
well to the other and the oscillation period T, namely, 

T= f d 1 = ^ A rcsin(^) (4.1) 

J -a V(Z) OJ \ H J 



f a + z , 1 2vr\/2 
T = 2 / dz 



l + ^Arcsinf^ 

z v{z) oj it \ p t 

= T + 2r (4.2) 

being To the period of the orbit of energy p for a particle with mass 2M in a single well of 
frequency u; one may notice that under the replacement p = M our, Fig. 6 describes as well 
the classical phase portrait of a particle in the ID-double well. We realize then that when 

6 



Vq is sufficiently low, the fluctuation essentially oscillates around the condensate along the 
classical orbit with period T ; in particular, when Vo equals //, the excitation requires twice 
a period to cross the double condensate, the classical orbit -pictures in medium dashes in 
Fig. 6 - being the separatrix between the two different oscillation regimes. As Vq grows 
up, a delay equal to 2 r appears in the oscillation time, reflecting the finite amount of time 
requested to switch back and forth between the two wells. 

The decrease of the eigenfrequencies can thus be unambiguously interpreted in terms 
of the enlarging of the classical period of the orbits due to the presence of a finite bar- 
rier. However, the most striking outcome of the TF computation is the discontinuity of 
the odd-state engenfrequencies as we approach the limiting regime n — Vq from above. In 
fact, contrasting this behavior with the strict Fourier spectrum, for which some results are 
indicated with symbols in Fig. 3, we realize that this discontinuity is an artifact of the TF 
approach, removable by an exact calculation. The explanation is as follows: let us project 
the wave equation onto an eigenstate, thus expressing the eigenvalues in the form 



05 = h ' ,h 



d5p h 



< I: 



[// - V ext (r)} (4.3) 



where the integration is to be performed within the borders of the TF condensate. If \x is 
lower than the barrier Vq, the condensate is split into two disjoint condensates and the region 
around z = is excluded from the integration domain. This exclusion is not significant for 
even parity states, since the density gradient vanishes at z = 0; however, this is not the case 
for the odd parity levels, whose exact Fourier spectrum, lying higher than the TF one (cf. 
Fig. 3) includes the contributions to the eigenfrequencies which stem from the whole spatial 
domain, in particular, the subbarrier region where the exact Gross-Pitaevskii wave function 
yields finite probability densities for the trapped bosons. 

This issue can be also easily interpreted taking into account that if one excites the 
condensate, generating a small density fluctuation with even parity, this excitation will 
oscillate in time within the limits of each well, basically unaffected by the presence of the 
other. By contrast, if the fluctuation is an odd parity one, a peak in one side of the condensate 
is imaged by a dip in the other; subsequent time evolution forces this pattern to oscillate 
from one side to the other, crossing the barrier. Now, since within the TF approximation 
both the condensate and its fluctuations are bounded to the region V ext < //, we recognize 
that when the barrier increases above the chemical potential value, any pertubation is forced 
to remain on its side with no posible interwell transition - thus the levels correspond to two 
independent Bose-Einstein condensates-. This is not what happens in a real double well 
described by the full GPE, where both condensates interpenetrate. 

The 3D situation gives us the opportunity to investigate transition density profiles for 
different selections of labels and for several possible manners of disecting the condensate. In 
Sec. 3 we have displayed some of the spatial characteristics of the first odd and even m=0 
and 1 modes (in Figs. 9 to 12), choosing the variable z to parametrize various equiprob- 
ability contours inside the condensate. The evolution of these patterns illustrates several 
interesting facts. We first realize that the spatial distribution of the perturbation evolves, 
with increasing well separation, so as to eliminate nodes of the eigenfunction, in order to ac- 
count for the level merging visible in Figs. 3 and 7. Secondly, we observe the transformation 
from a pattern localized in one harmonic potential to a pattern of the same nature, however 
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split into the two disjoint wells. In particular, the fourth portrait of each of these figures 
corresponds to an almost split condensate. We notice that in neither case the exact form of 
the twofold distribution after separation coincides with the original one, due to the fact that 
in the large Vq limit, the fluctuation is of the form (|2.11|) with the parity of Spo uncorrelated 
to that of the starting 5 p. Thus, should one duplicate the original pattern in each well after 
their definitive splitting, the energy would be higher. We observe that the m=0 modes dis- 
play finite density fluctuations in the outer edges of the condensate (z = ±(zq + Z)) whereas 
in z = 0, where the inner boundary surface grows up with increasing well separation, the 
density vanishes or not according to the parity being with odd or even; by contrast, in either 
case the m—1 modes intend to concentrate the fluctuations at values of z far from either the 
outer (i.e., z ~ ±(zo + Z)) or the inner [z ~ 0) condensate edges. The effect of the barrier 
upon even parity perturbations is to introduce a discontinuity in the slope of Sp, similar 
to that in the TF density itself. We also see that, as in the single- well problem, the inner 
contours, where the above density is larger, experience weaker density fluctuations than the 
condensate edge. 

V. SUMMARY 

In this work we have discussed the characteristics of the small amplitude collective ex- 
citations of a Bose condensed gas in harmonic wells where the confinement has undergone 
some splitting along one of the three directions. The excitations of the confined bosons has 
been investigated within the hydrodynamic limit of the GPE (TF limit), in other words, we 
have analyzed the small fluctuations of the condensate density in the limit of large number 
of particles. This enables us to determine the chemical potential as a function of the bar- 
rier height and the single-well parameters in the ID and 3D problems. We have considered 
extensively the ID potential consisting of two harmonic wells separated a distance 2a, since 
it contains essentially all the barrier effects that one may visualize later in the 3D situation. 

The general results can be summarized as follows. We find that the energy levels corre- 
sponding to the density fluctuations gradually merge into pairs as one increases the barrier 
height; within the TF approximation, as soon as the barrier overcomes the current value of 
the chemical potential the spectrum becomes completely degenerate. We also have shown 
that the behavior of the eigenspectrum as a function of the well separation can be inter- 
preted, in the ID case, by means of a classical picture, taking into consideration the way 
in which the finite barrier affects the period of the classical orbit. The appearance of de- 
generacy in the eigenspectrum, that images the well-known property of the Schrodinger 
eigenspectrum of a ID harmonic double potential, gives rise to various consequences upon 
the spatial pattern of the excitation near the surface of a 3D condensate when the latter has 
undergone splitting. In the 3D problem we have displayed typical excitations corresponding 
to levels with and without z-proyection of the angular momentum and have discussed the 
consequences of progressively splitting the condensate. We recognize the usual differences 
between fluctuations with zero and nonzero magnetic quantum number; while the first ones 
remain finite at the edges z = ±(z + %) (cf- Figs. 9 and 10), they vanish at those positions 
if m t^ 0. These features remain valid as one increases the barrier height and splits the 
condensate into two pieces. 

To summarize, we remark that the collective excitations of this bistable system can be 

8 



examined on almost identical grounds as in the case of a single-well one, with an additional 
parameter, namely the ratio between barrier height and chemical potential, which ranges 
from a single-well problem for negligibles barriers to the double degenerate regime where 
one exactly halves the condensate into two pieces, in the TF limit of the GPE for infinite 
separations. 
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FIGURE CAPTIONS 

FIG. 1. The chemical potential /i for the ID harmonic double well in the TF limit in 
units of the single well value /i as a function of the dimensionless well separation (in 
units of the size of the single- well condensate) . The horizontal and vertical thin lines 
respectively indicate the chemical potential and barrier height corresponding to full 
splitting into two disjoint single wells. 

FIG. 2. Three TF density profiles (in arbitrary units). Full lines indicate the value of 
the chemical potential for which the given contour is a condensate edge. The den- 
sity profiles have been moved up to the corresponding chemical potential for easier 
visualization. 



FIG. 3. Excitation spectrum (in units of u>) of the ID hydrodynamic equation (|2.3j ) as a 



function of the dimensionless parameter JVo/pL. The horizontal thin lines indicate the 
spectrum for the single well problem. Some eigenvales provided by the fit calculation 
are displayed with symbols. 

FIG. 4. The fast Fourier transform of the evolution of a perturbing density in the 1D-TF 



limit, for JVo/fi— 0, 0.4, 1.09 and 3.9 (from top to bottom). The intensity of the 
peaks is displayed in arbitrary units and Q is given in units of u. 

FIG. 5. Same as Fig. 1 for the 3D double well problem. 

FIG. 6. Three equiprobability contours (in arbitrary units) for a TF condensate, or equiv- 
alently, the condensate edges for the same chemical potentials considered in Fig. 2. 

FIG. 7. Same as Fig. 3 for the 3D hydrodynamic problem: (a) m — and (b) m = 1. 
Thin lines correspond to the single-well spectrum for the same parameters. 
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FIG. 8. Spectrum of the 3D problem provided by the fft calculation for wVb//i=0.14 as a 
function of Q (in units of to). The vertical lines show a few hydrodynamic eigenfre- 
quencies obtained through diagonalization of ( |2.2| ). The height of the peaks is set in 
arbitrary units. 

FIG. 9. The first odd-parity density fluctuation 5p[r(z), z] of the condensate as a function 
of z (in arbitrary units), for m=0. The four portraits show the evolution from a single 



well problem to the limiting case Vq ~ /i with the values JVof/j, = 0.01, 0.26, 0.43 
and 0.99. Each curve on each portrait corresponds to the contour values Po/p max = 0, 
1/4, 1/2, 3/4; the full lines denoting the condensate edge and ±(z + Z) points on the 
z-axis. 

FIG. 10. Same as Fig. 9 for the first even parity eigenstate. 

FIG. 11. Same as Fig. 9 for m=l. 

FIG. 12. Same as Fig. 10 for m=l. 
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